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Measurement of Interstitial 
Convective Heat Transfer 
and Frictional Drag for Flow 
Across Metal Foams 

Convective heat transfer and friction drag in a duct inserted with aluminum foams have 
been studied experimentally. The combined effects of foam porosity (e =^0. 7, 0.8, and 0. 95 ) 
and flow Reynolds number (I90O^Re^78OO) are examined. Frictional drags for flow 
across the aluminum foam are measured by pressure taps, while interstitial heat transfer 
coefficients in the aluminum foam are determined using a transient single-blow technique 
with a thermal non-equilibrium two-equation model. Solid material temperature distribu- 
tion is further measured for double check of the heat transfer results. To understand the 
frictional drag mechanisms, smoke-wire flow visualization is conducted in the aluminum- 
foam ducts. Results show that both the friction factor and the volumetric heat transfer 
coefficient increase with decreasing the foam porosity at a fixed Reynolds number. In 
addition, the aluminum foam of e =0.8 has the best thermal performance under the same 
pumping power constraint among the three aluminum foams investigated. Finally, empiri- 
cal correlations for pore Nusselt number are developed in terms of pore Reynolds number 
under various foam porosities. [DOI : 10.1115/1.141 6690] 

Keywords: Convection, Heat Transfer, Inverse. Non-Equilibrium, Porous Media 



Introduction 

Porous media have been widely used. in industrial applications 
such as electronic cooling, thermal energy absorber, geothermal 
system and many others. They can be in the form of packed beds, 
sintered materials, or foam materials. In the past, a large amount 
of research works for packed beds and sintered materials has been 
done [1-10], and extensive data and empirical correlations for 
convective heat transfer coefficient and friction factor are avail- 
able now, which cover a large range of particle diameters, Rey- 
nolds numbers and fluids. The heat transfer and friction in porous 
foam materials, however, has not been studied to the same extent 
as that of pack beds and sintered materials [11-14]. Thus, the 
related data and correlations are relatively sparse. It is inappropri- 
ate to use the data or correlations of packed beds and sintered 
materials to predict the convective heat transfer inside foam ma- 
terials because their internal structures are quite different [13]. In 
addition, the porosity of foam materials (e>0.7) is much larger 
than the porosity of the packed beds or sintered materials (e 
<0.476) [15]. Therefore, substantial data about the fluid flow and 
heat transfer characteristics inside foam materials are required to 
develop and optimize the design of such advanced devices. 

Two kinds of heat transfer coefficients are often employed to 
depict heat u-ansfer characteristics in porous media for their own 
application purpose, i.e., wall heat u-ansfer coefficient and inter- 
stitial heat u-ansfer coefficient. The first one always illusu-ates the 
global extent of heat transfer augmentation from the wall surface 
where the porous medium attaches. These surface heat u-ansfer 
data are useful in the design of heat sinks such as electronics 
cooling [8,10,11,16]. The second one describes the heat exchange 
between the fluid stream and the solid matrix of the porous me- 
dium. Technologic applications of these data include thermal en- 
ergy absorbers, advanced burners or incinerators, compact heat 
exchangers, etc. [7,12,13,17]. The present paper focuses on the 
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interstitial heat u-ansfer characteristics of the aluminum foams. 
Several relevant studies related to the interstitial heat transfer 
characteristics in porous media are briefly reviewed below. 
Golombok et al. [17] employed the heat regenerator technique to 
measure the heat transfer coefficient between the flowing gas and 
the metal fiber of a burner. The porosity of die porous fibrous 
material was fixed at e=0.8. Results showed that the heat u-ansfer 
coefficient for laminar flows increased rapidly with increasing gas 
flow. Vafai and Sozen [6] analyzed the thermal-fluid characteris- 
tics for flow through a packed bed of spherical particles of a fixed 
porosity e=0.39. They used the fluid-to-solid heal u-ansfer coeffi- 
cients from an empirical correlation established by Gamson et ai. 
[18]. Results showed that the local thermal equilibrium condition 
was very sensitive to the particle Reynolds number and the Darcy 
number. Younis and Viskanta [12] utilized a single-blow method 
to determine the volumeuic heat u-ansfer coefficients between the 
air stream and ceramic foam under a small range of the porosity 
(0.83 ^€^0.87). It was found that the volumetric heat u-ansfer 
coefficients were significantly different from previous data for 
packed beds and sintered metals. Ichimiya [13] evaluated the con- 
vective heat u-ansfer between the fluid and the solid material of an 
aluminum-ceramic foam (e=0.87) by comparing the measured 
and predicted Nusselt number on the heated wall of a porous 
channel. In his prediction, the analysis was performed under the 
assumption that porous media was consuoicted by small spherical 
particles. 

It is evident that the above studies related to the interstitial heat 
u-ansfer in porous media were performed under a fixed porosity or 
a small range of porosity. Suictly speaking, the results concerning 
the effect of porosity on the interstitial heat u-ansfer and friction of 
metal foam have not been reported yet. Therefore, the objecuve of 
this paper is to investigate the combined effects of porosity (e 
=0.7, 0.8 and 0.95) of the aluminum foam and the flow Reynolds 
number (1900^Re<7800) on interstitial heat transfer and friction 
characteristics in an aluminum-foam channel. Several issues are 
discussed in the present paper. First, experiments are conducted to 
examine the effect of porosity on the pressure-drop for flow across 
the aluminum foam, which has not been studied before. To quali- 
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tatively understand the frictional drag mechanism in the channel 
inserted with different-porosity aluminum foams, flow visualiza- 
tion using a smoke-wire technique is undertaken. Then, the volu- 
metric heat transfer coefficient between the gas and the aluminum 
foam is obtained by using a u-ansient single-blow technique. It is 
hard to directly measure the interstitial heat transfer coefficient of 
the porous medium because of the difficulty in the measurement 
of the temperatiu-e difference between the gas and solid phases. 
Therefore, an inverse method of transient single-blow technique is 
used to determine the convective heat u-ansfer coefficients be- 
tween the fluid and solid phases of the aluminum foams. Subse- 
quently, by quoting the friction and heat transfer data, a compari- 
son of the thermal performance among the three aluminum foams 
is made under two constraints, i.e., constant pumping power and 
constant though flow rate. Finally, empirical correlations for pore 
Nusselt numbers are developed in terms of pore Reynolds number 
for the different-porosity aluminum foam, which may be helpful 
in the design of related devices such as compact heat exchanger 
and thermal energy absorber. In addition, the interstitial heat u-ans- 
fer and friction data presented can provide a reference of 
computational-fluid-dynamic-based studies reladng the heat u-ans- 
fer in the foam materials. 

Experimental Program 

Apparatus. An experimental setup for the measiu-ement of 
interstitial heat U-ansfer characteristic in the aluminum foam is 
shown in Fig. 1. It consists of an airflow circuit, the test section, 
and insu-uments. Air from a 2 hp compressor first enters into a 
surge tank for the supply of airflow in the circuit. The compressed 
air released from the surge tank passes through two filters and a 
refrigeration dehydrator to remove oil, water and particulate ma- 
terial. It then flows into a setding chamber to reduce the possible 
noise and to attain a uniform axial velocity profile. Before enter- 
ing the test section, airflow is heated by passing through a spe- 
cially designed plate heater, which is formed by 30-Atni thick 
stainless steel foil wound around in a bakelite frame. The hot air 
then traverses the aluminum-foam test section, and finally exits 
from the test section to the outside of the building via an exhaust 
system. A digital flow meter is situated at the downsu-eam of the 
test section to measure the volumedic flow rate through the test 
section. The test channel as shown in Fig. 2 made of Plexiglas® 
plates has a rectangular cross section of 60x25.4 mm^. The alu- 
minum foam of size 60x25.4x60 mm^ is cut in the way that the 
sides are flush with the channel walls for a good contact. The 
aluminum foam is fitted into the channel, and subsequenUy the 
entire test section is bolted tightly from outside. Two pressure taps 
located at the inlet and the outlet of the test section measure the 
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Fig. 1 Sketch of the experimental setup 
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Fig. 2 Schematic drawing of the test section 



static pressure differences across the test sample. They are con- 
nected to a micro-differential transmitter to display the pressure 
signals. At both the upstream and the downstream of the test sec- 
tion, five diermocouple rakes (each with three beads) are used to 
measure the averaged airflow temperature. The transient responses 
of temperature signals are transferred and subsequently recorded 
by a real-time hybrid recorder (YOKOGAWA DAIOO) with a 
sample rate of 2 Hz. Since it is difficult to measure direcUy the 
solid-phase temperature of the porous matrix, twelve segments of 
copper plate of size 60X6X2 mm^ are fixed on the bottom surface 
of the foam metal for the measurement of wall temperature (Fig. 
2), which represents the solid-phase temperature of the porous 
maunx at the corresponding axial station. Using high-conductive 
epoxy at each interface between the aluminum fiber and copper 
plate minimizes the contact resistance. In addition, the copper 
segments are spaced by wood su^ips to minimize axial conduction 
between them. Twelve copper-constantan thermocouples are dis- 
tributed along the span-center of each segment, and their junction 
beads are carefully embedded into the copper segments. In addi- 
tion, to check the spanwise distribution of die wall temperature, 
the lateral wall temperature of the first, sixth, and twelfth seg- 
ments are measured by three thermocouples located at 20 mm 
away from the span-center of the segment. Results show that the 
wall temperature difference between the span-center and lateral 
thermocouples is less than 0.2®C. To ensure that the wall tempera- 
ture measured above can represent the corresponding solid-phase 
temperature in the porous matrix, an additional measurement of 
solid-phase temperature at the real face (x-L) of the foam ma- 
terial are conducted by infrared thermography. The infrared cam- 
era is focused on the exit plane of the porous medium from the 
channel outlet to view the solid-phase temperature disuibution. 
Results show that the solid-phase temperature disu-ibution is 
rather uniform, and the averaged solid-phase temperature compare 
well with the wall temperature with miaximum deviauon less than 
2°C. 

Porous Medium. Figure 3 is a photo showing the su-uctiue of 
aluminum foam (Duocel®, ERG Inc.) investigated, which is 
manufactured by directional solidification of metal from a super- 
heated liquid state in an environment of overpressure and high 
vacuum. It has a reticulated structure of open, duodecahedronal- 
shaped cells connected by continuous, solid metal ligaments. The 
matrix of cells and ligaments is completely repeatable, regular, 
and uniform throughout the entirety of die material. The solid 
ligaments are made of aluminum alloy 6101-T6 that has a thermal 
conductivity of 218 W/m-K. Table I summarizes the properties of 
the three samples of the aluminum foams that are used in the 
experimental study. The effective thermal conductivity of the po- 
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Fig. 3 Photograph of the aluminum foam matrix 



rous matrix is expressed by a function of geometry and solid 
conductivity, i.e., k* = (l -e)k,. The volume-based porosity (e) 
provided by the manufacturer is deiermined using the weight of 
the sample widi the density of the aluminum. It is rechecked by 
determining the void-to-total volume ratio of the sample, in which 
the sample total volume is measured by its geometry and the pore 
volume is determined by filling water inside the sample. It should 
be noted that, due to the manufacturing limitation, the aluminum 
foam produced by ERG has a small range of porosity about 
0.88^6^0.97. In this study, the low-porosity samples of e=0.7 
and 0.8 are obtained by slowly compressing the aluminum foam 
of 6=0.95 (pore density 10 PPI) widi identical compression ratio 
in all directions [19,20]. Measurements of pressure drop/thermal 
conductivity by rotating the porous matrix relative to flow/ 
temperature gradient direction indicate that the mau-ix is 
hydraulically/thermally isotropic. In addition, visual inspections 
do not reveal any non-uniformity in the su^cture of matrix. As for 
the permeability and inertia coefficient of the three aluminum 
foams, they are deiermined by the fully developed pressure-drop 
across the porous sample together with the mean velocity in the 
porous channel under an isothermal condition. A rearranged ver- 
sion of modified Darcy equation [11] expresses the relationship 
between these two constants. 



PfU 



(1) 



Figure 4 shows the relationship between (^P/L)/ipLu) and 
pfu/jLL for three different porosities investigated. The su-aight lines 
are curve-fit results of these measured data with standard devia- 
tion less than 1 1 percent. The permeability K and inertial coeffi- 
cient F are determined from the intercept together with the slope 
of the line [11]. According to the estimation mediod of Antohe 
et al. [19], die uncertainties of AT and F are less than 6 percent and 
9 percent, respectively. Note that the length of the present porous 
matrix should be long enough such that die value of ^P/L can 
well represent the fully developed pressure gradient dP/dx in Eq. 
(1). To this end, experiments of longer porous matrices of L 
= 90 mm and 120 mm, are conducted to check the enu-ance and 
exit effects on die pressure gradient- Results show that in the 
range of flow velocity investigated the value of AP/L for L 
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=60 mm is about 5.6 percent and 5.3 percent higher than those of 
L = 90 mm and 120 mm, respectively, which represents the typical 
entrance and exit effects of the present porous matrix. 

Theoretical Simulation 

Heat transfer in a porous channel involves a complex thermal 
transport mechanism between the air and surface of solid matrix. 
It is difficult to find the exact solutions for the flow and tempera- 
ture fields. The U-ansient single-blow technique utilized in the 
present study offers a hot fluid to heat up the porous medium. The 
physical configuration considered is a rectangular channel filled 
with the aluminum foam. High temperature airflow passes through 
and transports energy to the porous medium. The heat transfer 
characteristic of die porous medium can be evaluated by the 
transient response of die ouUei air temperature, and double re- 
checked by die wall temperature disuibution. To establish a math- 
emadcal model, die following assumptions and simplifications are 
employed. 

1 The inlet fluid temperature is uniform. Temperature varia- 
tions of the solid and fluid in the porous channel depend only 
on die position along the flow direction and time; therefore, it 
is a one-dimensional u-ansient problem. 

2 The medium is isotropic and has a uniform porosity. 

3 The fluid flow is steady and incompressible; its thermal 
physical properties are constant. 

4 There is no local thermal equilibrium between air and solid. 

5 The fluid velocity in the porous channel is uniform. 

The governing equauons are developed by application of the local 
volume-average technique. This averaging process may obscure 
local pore phenomena that contribute to the global u-ansport. En- 
ergy equauons for die solid and fluid phases are as follows: 



Table 1 Physicai properties of the aluminum foam 



Porosity, e 


0.70 


0.80 


0.95 


Average Ligament Diameter, La (mm) 


0.36 


0.36 


0.36 


Thermal Conducdvity, ik/ (W/m-K) 


65.4 


43.6 


10.9 


Mean Cell Size, d (mm) 


1.84 


1.92 


2.03 


Permeability, JCxlO*(m^) 


0.88 


2 


7.5 


Inertial Coefficient, F 


0.11 


0.068 


0.065 
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Table 2 Coefficients under various fiow parameters 
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0.7 


1900 


8.2 


4.0 


2800 


8.8 


4.7 


3800 


8.1 


3.2 


6100 


8.8 


2.1 


7800 


9.5 


1.8 


0.8 


1900 


4.7 


0.1 


2800 


6.0 


0.8 


3800 


6.5 


0.1 


6100 


7.7 


O.l 


7800 


7.6 


0.1 


0.9 


1900 


11.1 


0.2 


2800 


13.9 


1.6 


3800 


12.1 


1.3 


6100 
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(I - c )( 1 + ?, )(pcp),-5ji = /.„(r/- T,) + 



(2) 



the governing equations can be nondimensionalized as; 



I dTf U dTf\ ^ ^l^^Tf\ 

In Eq. (2) fi is a constant to modify the time delay due to heat 
capacity and heat loss of the copper plate on the porous medium. 

and /:* are the equivalent conductivity of the solid and fluid 
phases, respectively, and are represented as ekf and (l-e)A:^. 
Moreover, k* is the transverse thermal dispersion conductivity 
due the fluctuating tortuous air motion in the solid matrix. Accord- 
ing 10 the simple model reported by Hunt and Tien [11], die dis- 
persion conductivity is expressed as k* -0.025(pCp)fyjKU. 
The appropriate initial and boundary conditions are 



T,{0^) = Tf(0^) = Tr 



(4) 



( 1 - e)^.-^ (^o) = ( I + (2)hMTi{t) - r,(/.o)] (5) 



Tf{t,0) = TM 



(6) 
(7) 
(8) 



^2 is another constant to modify the inlet boundary condidon due 
to the heat u-ansfer enhancement by the random structure of the 

porous medium. These weighting factors could be determined by 
trial and error via comparing the measured and predicted wall 
temperature distributions. Table 2 summarizes the values of d 
and f 2 selected for different flow and parameters in the present 
study, r, is die inlet air temperature and is determined from the 
measurement. By introducing the following variables: 

^_ T-Tr t 
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C2 d^e, 
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(10) 



The dimensionless forms of the initial and boundary conditions 
are as follows: 



<l,(0^)=fl/(OJf) = 0 
(1 + WB 



1-e 



[e,-S,(r,0)]=O 



d9. 
ddf 



(11) 

(12) 

(13) 
(14) 

(15) 



Using a guessed value of volumetric heat transfer coefficient (/i^) 
together with the prescribed values of geometric properties, physi- 
cal properties, measured inlet air temperatures, a volumetric flow 
rate, and proper choices of empirical constants of and (2^ the 
solid and fluid phase temperatures, i.e., 0^ and Of, can be solved 
from Eqs. (9) and (10). Then, the value of /i„ is modified until the 
accuracy between predicted and measured exit air temperatures is 
satisfactory (typically with deviation less than l.OxlO"^). 
The flowchart of the transient single-blow technique is showed in 
Fig. 5. 

The methodology used to solve the above equations associated 
with the boundary and initial condidons is based on the control 
volume, finite-difference technique [21]. The u-ansient finite- 
difference form of the energy equations is solved expliciUy by a 
line-by-line iterative method. Stepping forward in time and retain- 
ing a converged solution at each time step handle the time depen- 
dence. The sensitivity studies of node number and time increment 
are performed to check the stability and accuracy of die numerical 
scheme. Five sets of node numbers, i.e., 40, 60, 100, 150, and 200, 
are tested to ensure that the results of die present study are grid 
independent, A comparison of Of between 100 nodes and 150 
nodes reveals accuracy up to three significant figures, suggesung 
that 100 nodes (AX = 0.01) are sufficient to obtain grid indepen- 
dent results. In addition, the effect of time step on the results is 
examined by testing several time steps of Ar=0.05, O.l, 0.2, and 
0.5. The results indicate that there is no significant difference be- 
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Fig. 5 Flow chart of the transient single-blow technique 



iween the results using time increments of 0.05 -and 0.1. Hence, 
the entire calculations are carried out with AX=O.Ol and 

At =0.05. 

Data Analysis and Uncertainty 

Friction Factors. The friction factor of the present foam 
channel is calculated from the pressure drop across die test chan- 
nel and the bulk velocity of the air and expressed as 

(16) 

Note that based on the heating level of this study, it is experimen- 
tally determined that the friction factor with heating is only about 
2 percent higher than that without heating. Therefore, for time 
saving, the friction factor is based on the isothermal conditions 
without the tedious heating process. The individual uncertainties 
for AP and C/ are 5.5 percent and 3.0 percent, respectively. The 
maximum uncertainty of Cf is estimated to be less than 7.3 per- 
cent for Reynolds number greater than 1900 by the uncertainty 
estimation method of Kline and McClintock [22]. 

Interstitial Heat Transfer Coefficient. The errors of heat 
transfer coefficient are produced mainly due to the deviation of 
temperature measurement and the root-mean-square of the differ- 
ence between the values of experiment and computation. The un- 
certainty in temperature measurement is about 3.3 percent. The 
maximum uncertainty in die root-mean-square of the temperature 
difference is less 14.0 percent. The combination of root-sum- 
square is employed to calculate die overall uncertainty of die 
volumeuic heat transfer coefficient. As a result, die overall uncer- 
tainty of /i„ is 14.0 percent. 



Results and Discussion 

Friction Factor. The effect of the Reynolds number (Re) on 
Uie friction factor (C/) for die flow across die channel inserted 
wiUi die aluminum foam is shown in Fig. 6. Three aluminum 
foams of different porosities are tested, i.e., e=0.7, 0.8, and 0.95. 
The data points are die exact experiment and die lines passing 
through diese points are least-square results. It is seen dial the 
friction factor decreases widi increasing Reynolds number for all 
porosities investigated. The aluminum foam of higher porosity has 
a lower friction factor due to the less channel blockage togedier 
widi a smaller surface drag. To qualitatively understand die drag 
of flow across different-porosity foams, a comparison of die flow 
patterns around different-porosity aluminum foams is shown in 
Fig. 7, which is measured by a smoke-wire technique [23]. The 
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Fig. 6 Friction factor in the aluminum foam channel 



smoke wire is placed vertically along die duct center on the axial 
station ahead of aluminum foam. Instead of die diick-sample (60 
mm) employed in die heat u-ansfer experiment, a slender alumi- 
num foam (10 mm) is used herein for easier permeable. The air- 
flow is from left to right. Figure 6(a) shows the flow over a 
smooth plane channel without porous medium and serves as a 
reference. Figures 6(b) and (c) are die results of e=0.7 and 0.95, 
respectively. For e=0.7, it is seen dial some smokes stay in front 
of die aluminum foam, meaning that the flow cannot easily pen- 
etrate die aluminum foam due to a large blockage. In conu-ast, the 
flow for £=0.95 seems to pass smoodily through die aluminum 
foam. The above flow-structure comparison has reasonably illus- 
trated die effect of porosity on die friction factor showing in Fig. 
6. Anodier interesting observation by comparing Figs. 1(b) and (c) 
is that die flow suiictures behind die aluminum foam are different. 
A large amount of eddies is clearly observed behind the rear face 
of the aluminum foam for e=0.7 and. however, the flow after die 
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Fig. 7 Comparison of flow structure around the aluminum 
foams of different porosities 
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Table 3 Coefficients of friction factor correlations of the porous channel 



e 






ai 


Deviations 


Eq. 


0.95 


141 


-0.397 


4.2% 


(17-1) 


0.80 


160 


-0.562 


5.5% 


(17-2) 


0.70 


327 


-0.546 


3.4% 


(17-3) 



high-porosity aluminum foam (e=0.95) is relatively laminar. 
Therefore, we may extrapolate this fact to that the flow inside the 
foam material of e=0.7 would be more turbulent than that of 
e=0.95, which largely affects the heat transfer characteristics in- 
side the porous medium and will be shown later. The data pre- 
sented in Fig. 5 could be correlated as the following equations: 

Q=aiRe^^ (17) 

where the coefficients ai and a 2 for various foam porosities are 
listed in Table 3. The maximum deviation between the equation 
above and the experimental data is only 5.5 percent. 

Heat Transfer Characteristics. TVpical examples of the his- 
tories of inlet and outlet air temperatures during the present ex- 
periment for the single-blow technique are shown in Fig. 8 under 
various e and Re. The solid circles in each figure are the measured 
inlet air temperatures, while the solid lines passing through these 
circles are the curve-fitted results (by the 5th order polynomials) 
that serve as the inlet condition to the present theoretic model 
equations. The solid u-iangles are the measured outlet air tempera- 
tures, while the dashed lines are the calculating results from the 
present theoretic model using the volumetric heat transfer coeffi- 
cient listed in each figure. In general, the theoretically simulated 
outlet air temperatures are in good agreement with the experimen- 
tal ones, indicating the appropriate volumetric heat u-ansfer coef- 
ficient selected in each case. It should be noted that in all tests the 
heat dissipated from the electric heater is controlled to be constant 
by adjusting the current through the resistance of the heater. Al- 
though the range of abscissa (time) displayed in each figure may 
be different, the histories of the inlet air temperatures under the 
same mass flow rate (Re) are essentially the same. As shown in 
Figs. 8(fl)-(c) for Re= 1900, with the same history of the inlet air 
temperature, the smaller the foam porosity is, the slower is the 
incremental rate of the outlet air temperature. For example, at time 
step of r = 100, the outlet air temperatures are So-O.Oll, 0.016, 
and 0.008, respectively, for e=0.95, 0.08, and 0.7. This means 
that the aluminum foam of a small porosity has a large heat ca- 
pacity that reduces the incremental rate of the outlet air tempera- 
ture. There is a significant heat exchange between the air and solid 
matrix for the small- porosity aluminum foam. Attention is now 
turned to the effect of the Reynolds number showed in Fig. 8(c), 
in which the difference between the inlet and outiet air tempera- 
tures is roughly the same for Re=1900 and 6100. It has clearly 
showed that the higher Reynolds number (Re=6100) requires 
more time to reach the same temperature difference than the lower 
Reynolds number (Re=l900). This is because a higher Reynolds 
number has^a higher heat transfer coefficient inside the porous 
medium that removes more heat from the air to the porous me- 
dium, and thus requires more time to increase the difference of the 
inlet and outlet air temperatures. Figure 9 presents the temperature 
disuibution of the solid material (the aluminum foam) along the 
streamwise direction. It is found that the theoretic curves agree 
fairly well with the experimental data. This confirms again the 
correctness of the predicted heat transfer coefficient. 

There is a broad variety of constraints that may be utilized in 
evaluating the thermal performance of the channel with porous 
medium [24,25]. In this paper, thermal performance comparisons 
among the three aluminum-foam channels are made for two sets 
of constraints that are widely employed, namely, the same through 
flow rate and the same pumping power. The comparison of heat 



u-ansfer coefficients for the same airflow rate consu-aint is made in 
Fig. 10, which shows the Nusselt number as a function of the 
Reynolds number for three different aluminum foams. The corre- 
sponding volumeuic heat transfer coefficient (h^) is also plotted 
on the right-hand side of ordinate of this figure. As shown in, Fig. 
10, the Nusselt number increases with increasing the Reynolds 
number. For a fixed Reynolds number, the heat u-ansfer coefficient 
increases with decreasing the porosity due to an increase in local 
fluid velocity and, partly, the enhancement of turbulent u-anspori 
inside the foam (Fig. 7). Attention is now turned to die same 
pumping power constraint. The pumping power required to drive 
a fixed volumetric flow rate ( V) in the foam duct can be written as 
diP V. As shown in Fig. 6, the pressure drop is higher for the 
foam channel of a lower porosity. Therefore, the use of low- 
porosity foam instead of high-porosity foam should reduce the 
through flow rate to keep the fluid pumping power constant in a 
porous channel, which lead to a reduction in the heat transfer 
coefficient. To implement this case. Fig. 11 shows the Nusselt 
number as a function of the nondimensional pumping power Cf 
• Re^. From this figure, it is seen that the volumeuic heat transfer 
coefficient increases with increasing the nondimensional pumping 
power for all aluminum-foam channels. Further inspection of this 
figure reveals that the aluminum-foam channel of e=0.8 has the 
best thermal performance among the three aluminum-foam chan- 
nels. This is because the aluminum-foam channel for e=0.8 is 
accompanied by a moderate pressure-drop (Fig. 6) and produces a 
sufficiently high heat transfer, and thus yields the highest thermal 
performance. 

Comparison With Previous Results. A comparison of the 
present and previous interstitial heat u-ansfer results is made in 
Fig. 12. The volumeu-ic heat u-ansfer coefficient is expressed in 
dimensionless form in terms of the pore Nusselt number ( Nuj) , 
which is plotted against the pore Reynolds number (Re^) in Fig. 
12. In general, the mean pore diameter appears to be a reasonable 
choice for the characteristic length in defining the Nusselt number 
and Reynolds number in the foam materials [12]. In this figure, 
the present data are presented by open symbols, and the su-aight 
lines passing through these symbols are least-square fits, and can 
be expressed as 

m^=byRe^/. (18) 

The coefficients and bi for various porosities and the maxi- 
mum deviation between the equation above and the experimental 
data are listed in Table 4. Note that the data of pore Nusselt 
numbers for e=0.7 and 0.8 gather closely, and can be further 
correlated as a single equation with a reasonable deviation. As for 
the previous interstitial heat Uransfer results, the data band repre- 
sents the results obtained from ceramic foams of 0.83^c^0.87 by 
using a single-blow technique [12], which consists of three corre- 
lations for various pore diameters d = 0J6, 0.94 and 1.52 mm. 
The Reynolds-number dependence (Re J) for these correlations 
ranges from rt = 0.7 to 0.96. In addition. The dashed line is the 
volumetric heat u-ansfer coefficient inside a ceramic-aluminum 
foam for e=0.87, which is obtained by comparing the measured 
and predicted Nusselt number on the heated wall of a porous 
channel [13]. A light irregularity of the curve at Rej=20 is be- 
cause of the inconsistency of results obtained for rf=2.0 and 4.2 
mm at the same Reynolds number. Furthermore, the centeriine is 
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an empirical correlation compiled from transient and steady state 
data of packed spherical particles. This correlation is unique and 
covers a broad range of porosity and Reynolds number [25]. 

It is seen from this figure that the Reynolds-number dependence 
(ReJ) of the present data = 0.6-0.83) is largely similar to that 



of Younis and Viskanta [12]. Quantitatively, the present data for 
e=0.7 and 0.8 can exu-apolate well to those in low pore Reynolds 
numbers obtained by Ichimiya [13] but are slightly lower than 
those of Younis and Viskanta [12]. This discrepancy in heat u-ans- 
fer coefficient may be aiuibuted to the different size (in su-eam- 
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wise direction) of foam used between Younis and Viskanta [12] 
and the present woric. In Younis and Viskanta [12], a shorter ce- 
ramic porous medium (12-14 mm) than the present one (60 mm) 
will cause a more significant entrance effect that increases the 
length-mean heat u-ansfer coefficient. Nevertheless, the above 
comparison between the present data with those by Younis and 
Viskanta [12] and Ichinoiya [13] is satisfactorily in agreement. It is 
further seen that the correlaiion proposed by Wakas and Kaguei 
[26] seems to be invalid for the data presented, especially in the 
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low Reynolds number range. We attribute this to the fact that the 
porous structures of the present foam material are different from 
those of the packed beds, which has a large influence on the con- 
vective heat transfer through affecting the shape of the flow pas- 
sages as wet! as the surface area per unit volume. 

Conclusions 

Interstitial heat transfer and friction characteristics for flow 
across a channel inserted with alimiinum foams have been studied 
experimentally. The combined effects of foam porosity {e=0.70, 
0.80, and 0.95) and flow Reynolds number (1900^Re^7900) are 
examined in detail. A transient single-blow technique has been 
applied to the determination of the interstitial heat transfer coef- 
ficient of the aluminum foam. The results of numerical simulation 
coincide with those of experimental measurement by using the 
thermal non-equilibrium energy equations together with appropri- 
ate thermal boundary conditions and proper choices of the empiri- 
cal constants d and (2- Main findings based on the data pre- 
sented are as follows. 

1 At a fixed Reynolds number, as the porosity of the aluminum 
foam increases from e=0.70 to 0.95, the friction factor decreases 
due to the decrease in the channel blockage. The friction factor 
decreases with an increase of Reynolds- number for all porosities 
investigated. 

2 Under the same through flow rate and history of the inlet air 
temperatiu-e, the small-porosity aluminum foam has a slow incre- 
mental rate of the outlet air temperatures, meaning that the alumi- 
num foam of a small porosity has a large heat capacity. 

3 At a fixed Reynolds number, the volumetric heat u-ansfer 
coefficient inside the aluminum foam increases with decreasing 
the porosity due to die increase in local fluid velocity and the 
enhancement of turbulence transport. In addition, the volumetric 
heat transfer coefficient increases widi increasing the Reynolds 
number under a fixed value of foam porosity. 

4 Among the three aluminum foams investigated, the foam 
channel of e=0,8 appears the best thermal performance under 
constant pumping power consu-aint because it produces a satisfac- 
torily high heat transfer and is accompanied by a moderate pres- 
sure drop. 

5 Empirical correlations of pore Nusselt number for three alu- 
minum foams under a specific length {t = 60mm) are developed 
in terms of pore Reynolds number for the first time, which can be 
helpful in the design of the related devices such as compact heat 
exchangers or thermal energy absorbers. Further efforts are re- 
quired to conclude a unique correlation that covers a wide range 
of the flow and geomeuic parameters such as the porosity of foam 
metal, Reynolds number, etc. 
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Nomenclature 

Bi = Biot number, h^L^/k, 

C| = constant C, = (pCp)^/[(pc^),(lH-f,)]; 

C2 = constant C2 = Mp<:p)//[*/pCp),(i+f,)] 

Cf = friction factor, Eq. 16 

Cp = isobaric specific heat, kJ/kg-K 

d = mean pore diameter, m 

F = inertia coefficient constant 
Ay = volumeuic heat transfer coefficient, W/m^-K 

K = permeability, m^ 

k = thermal conductivity, W/m-K 

L = test section length, m 
Nu = Nusselt number based on porous matrix length. 
h,L^/k^ 



Nurf = pore Nusselt number. h„d}ikf 
P = pressure, N/m^ 
Pr = Prandtl number {^iCp)flkf 
Re = Reynolds number, PfUU^ 
Re^ = pore Reynolds number. pjUdf{efj.) 
St = Stanton number, h^,L/[(pCp)fU] 
T = temperature, K 
T'finaj = steady state fluid temperature at test section inlet, K 
Tr = reference (room) temperature, K 
t = ume, s 

U - mean velocity in the rectangular channel, m/s 

u = pore velocity, m/s 

X = dimensionless axial length, x/L 

X = axial length, m 

e = porosity 

0 = dimensionless temperature (r-T^VCrgnai-r,) 
fjL = viscosity, kg/m-s 

p = density, kg/m^ 
T = dimensionless time, U t/L 
(1,(2 = constants 

Subscripts 

d = pore 

/ = fluid 

1 - inlet value 
s = solid 

Superscripts 

* = equivalent value 
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